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ABSTRACT 

We present new high resolution {R ~ 18, 000) near-infrared spectroscopic observa- 
tions of a sample of classical FU Orionis stars (FUors) and other young stars with FUor 
characteristics that are sources of Herbig-Haro flows. Spectra are presented for the 
region A = 2.203 - 2.236 fim which is rich in absorption lines sensitive to both effec- 
tive temperatures and surface gravities of stars. Both FUors and FUor-like stars show 
numerous broad and weak unidentified spectral features in this region. Spectra of the 
2.280 - 2.300 nm region are also presented, with the 2.2935 iim v=2-0 CO absorption 
bandhead being clearly the strongest feature seen in the spectra all FUors and Fuor-like 
stars. A cross-correlation analysis shows that FUor and FUor-like spectra in the 2.203 
- 2.236 nm region are not consistent with late-type dwarfs, giants, nor embedded proto- 
stars. The cross-correlations also show that the observed FUor-like Herbig-Haro energy 
sources have spectra that are substantively similar to those of FUors. Both object 
groups also have similar near- infrared colors. The large line widths and double-peaked 
nature of the spectra of the FUor-like stars are consistent with the established accretion 
disk model for FUors, also consistent with their near-infrared colors. It appears that 
young stars with FUor-like characteristics may be more common than projected from 
the relatively few known classical FUors. 
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Introduction 



FU Orionis stars (FUors) are rare. Classical FUors have been identified by large (~ 5 mag) 
increases in brightness and luminosity followed by fading over decades, and only a few young 
stars have been co nfirmed as such, iiotably FU Ori, V1057 Cyg, V1515 Cyg, V17 35 Cyg (jHerbig 
1973), V346 Nor (|Graham fc Frogei llQsd ^. and V733 Cep feeipurth et al.1 120071 ^. Addit ionally 
while in outburst, they exhibit optical and near-IR spectra similar to FU Orionis itself (jHerbig 



19771 : iHartmann &: Kenvon 



1996 



Herbig et alj 120031 ) . The absorption line widths (interpreted as 



rotational velocities) and derived spectral types of FUors change with wavelength, varying from 
F or G in the visible to M type in the near-IR. Their spectra also indicate supergiant or gi- 
ant surface gravities. These characteristics have been modeled as arising from young stars with 
massive accretion disks, th e latter dominating the spectra and luminosities of these objects (see 
Hartmann &: KenvonI 119961 ). The accretion disk model has been quite successful in explaining the 
basic visible-to-IR s pectral features and energy distributions of FUors (jHartmann &: Kenvonlll985l : 
Kenvon et al.lll988l: iHartmann et al.ll2004l : iGreen et al.ll2006l ). although in detail there are still dis- 
crepancies (jHerbig et al.ll2003l ). The only published (and modeled) high resolution near-IR spectra 
of FUors have been in the vicinity of the first vibrational overtone t; = — 2 of CO (A > 2.294/um), 
providing temperature and rotation information. New observations at other near-IR wavelengths 
may aid in determining the range of effective temperatures and surface gravities over which FUor 
spectral features arise in the near-IR, providing more constraints on applicable physical models. 

The general symmetry and relatively even spacing often seen in Herbig-Haro object knots 
strongly sugg est that they are somehow linked to episodic outbursts from their parent star plus 
disk systems (|Dopitalll978l : lReiDurthlll989l'l. This suggests th at the sources of HH objects and FUor 
outbursts may possibly be related. iReipurth AspinI (jl997l ) found that five young stars associated 
with HH flows (termed Herbig-Haro energy sources or HHENs) had low resolution K-band spectra 
that are very similar to those of FUors. If episodic outbursts do generally indicate FUor activity, 
then there may be many more young stars with FUor characteristics than the number predicted 
by extrapolating numbers from the few known classical FUors. In the following discussions, we use 
the term FUor-like star to indicate objects th at have spectral simi larities to the classical FUors, 
but for which no eruption has been witnessed (jReipurth et al.ll2002l ). 



How similar are classical FUors and FUor-like stars? Both are young stars that are likely 
undergoing episodes of high accretion. However, their luminosities may be powered by different 
physical processes or they may be at different phases of their outburst and decay cycles. It is also 
important to understand how similar these objects are to embedded (Class I) protostars accreting 
at fairly high rates (e.g., M ~ 10~^Mq yr~^) and to investigate whether there are any spectral 



^Much of the data presented herein were obtained at the W.M. Keck Observatory from telescope time allocated to 
the National Aeronautics and Space Administration through the agency's scientific partnership with the California 
Institute of Technology and the University of California. The Observatory was made possible by the generous financial 
support of the W.M. Keck Foundation. 
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similarities between protostars, FUors, and FUor-like objects. Also, do any of these objects show 
spectral characteristics similar to other young eruptive variables such as EX Lupi-type (EXors) 
young stars? EXors undergo multiple short duration optical outbursts of 1 - 5 mag but have K or 
M dwarf optical sp ectra dominated by emission lines, much different from FUor spectra (e.g., see 



Herbig et al.ll200lh 



We believe that such fundamental questions can be addressed with new high resolution spectra 
Such spectra must cover near-IR wavelengths since many FUor-like stars are highly extinguished 



and c annot be observed in visible light. The low resolution near-IR spectra of iReipurth &: Aspin 



(jl997l ) did show strong similarities between FUors and FUor-like objects, but those spectra were 



also consistent with giant stars. One object with such a nea r-IR spectrum in. close proximity to a 



young star was recently found to be a background giant star (jAspin &: Greendl2007l ). Linking other 
young stars to FUors will improve our knowledge of the number and types of young stars with 
high accretion rates, and detailed modeling may reveal more information on physical structures 
and accretion mechanisms or processes. 

We have conducted a new near-IR spectroscopic study of FUors and FUor-like stars covering 
a relatively broad range of wavelengths near 2 ^m, including spectral features sensitive to both 
effective temperature and surface gravity. We present these new data in §2. The properties of 
the individual FUor-like stars are discussed in §3. In §4 we analyze the similarity of the spectra 
of FUor-like stars to those of i) classical FUors, ii) late-type spectral standards, and in) Class I 
protostars. We discuss the likely physical similarities of these stars and the possible origins of their 
attributes in S5. 



2. Observations and Data Reduction 

High resolution near-IR spectra of several classical FUors and FUor-like stars were acquired 
with the NASA IRTF, Keck II, and GEMINI-South telescopes. The observation dates, total inte- 
gration times, signal-to-noise, equipment used, and coordinates of all observed objects are given in 
Table m 



2.1. Object Sample 



The observed FUors and FUor-like objects are listed in Table [H The FU or-like Herbig-Haro 
source s V883 Ori, HH 354 IRS, HH 381 IRS, and L1551 IRS 5 were observed bv beipurth Aspin 
(jl997l ) who found that their low resolution (R ~ 420) A'-band spectra were similar to those of FUors. 
Parsamian 21 has not to our knowledge been observed previously with near-infr ared spectroscopy. 
The spectra of L1551 IRS 5 were previously published in lDoppmann et al.l (|2005l ). New and unique 
observations were made of the FUor V1057 Cyg in the 2.2935 ^m C O bandhead region, a nd we 
use its previously published spectrum in the 2.2075 ^m Na line region ([Greene &: Ladalll997l ). New 
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spectra of FU Ori itself were also acquired. 



2.2. IRTF Observations 

Near-IR spectra of V1057 Cyg were acquired on UT 1999 August 30 January with the 3.0 m 
NASA Infrared Telescope Facili ty on Mauna Kea, Hawaii, using the CSH ELL facility single-order 



cryogenic echelle spectrograph (ITokunaga et al.lll99d : iGreene et al.lll993l ). Spectra were acquired 



with a I'.'O (5-pixel) wide slit on the dates indicated in Table [TJ This provided a spectroscopic 
resolution R = \/6X = 21,000 (14 km s~^). The spectrograph was fitted with a 256 x 256 pixel 
InSb detector array, and custom circular variable filters (CVFs) manufactured by Optical Coating 
Laboratories Incorporated were used for order sorting. These filters successfully eliminated the 
significant interference fringing normally produced in CSHELL and other echelle spectrographs 
which use CVFs for order sorting. The plate scale was 0'.'20 pixel" ^ along the 30" long slit (oriented 
east - west on the sky), and all spectra were acquired at a central wavelength setting of 2.29353 fim 
corresponding to the v=2-0 CO bandhead. Each exposure had a spectral range AA ~ A/400 {Av ~ 
700 km s"^). 

Data were acquired in pairs of exposures of up to 180 s duration each, with the telescope 
nodded 10" east or west between exposures so that object spectra were acquired in all exposures. 
The total integration time on V1057 Cyg was 13.0 minutes. The AlV star HR 8585 was observed at 
nearly identical airmass for telluric corrections. Spectra of the internal CSHELL continuum lamp 
were taken for flat fields, and exposures of the internal CSHELL Ar and Kr lamps were used for 
wavelength calibrations. 



2.3. Keck Observations 

Spectra of FU Ori and all five FUor-like stars were acquired on UT 2001 July 7-8, 2001 
November 6, and 2007 March 6. These data were obtained with the 10-m Keck II telescope 
on Mauna Kea, Haw aii, using the NIRSPEC multi-order cryogenic echelle facility spectrograph 



([McLean et al.lll998l ). Spectra were acquired with a 0"58 (4-pixel) wide slit, providing spectroscopic 
resolution R = X/5X = 18,000 (16.7 km s~^). The plate scale was 0'.'20 pixer^ along the 12" slit 
length, and the seeing was typically 0'.'5-0'.'6. The NIRSPEC gratings were oriented to allow orders 
containing the 2.1066 fim Mg and 2.1099 /xm Al lines, the 2.1661 //m HI Br 7 line, the 2.206 and 
2.209 fim Na lines, and the 2.2935 /xm CO bandhead regions to fall onto the instrument's 1024 x 
1024 pixel InSb detector array. The NIRSPEC-7 blocking filter was used to image these orders on 
the detector. NIRSPEC was configured to acquire simultaneously multiple cross-dispersed echelle 
orders 32-36 (2.08-2.37 /im, non-continuous) for all objects. Some objects were also observed in 
adjacent orders, resulting in a somewhat increased spectral range at longer or shorter wavelengths. 
Each order had an observed spectral range AA ~ A/67 {Av ~ 4450 km s~^). 
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The slit was held physically stationary during the exposures and thus rotated on the sky as 
the non-equatorially-mounted telescope tracked when observing. Data were acquired in pairs of 
exposures of durations from 5-500 s each, with the telescope nodded 6" along the slit between frames 
so that object spectra were acquired in all exposures. Early-type (B9-A2) dwarfs were observed 
for telluric correction of the FUor and FUor-like spectra. The telescope was automatically guided 
with frequent images from the NIRSPEC internal "SCAM" IR camera during all exposures of more 
than several seconds duration. Spectra of the internal NIRSPEC continuum lamp were taken for 
flat fields, and exposures of the Ar, Ne, Kr, and Xe lamps were used for wavelength calibrations. 



2.4. GEMINI Observations 



Spectra of FU Ori were acquired on UT 2006 April 03 with the Phoenix near-IR spectrograph 
(jHinkle et al.ll2002l ) on the 8-m GEMINI-South telescope on Cerro Pachon, Chile. Spectra were 
acquired with a 0'.'35 (4-pixel) wide slit, providing spectroscopic resolution R = X/5X = 40,000 
(7.5 km s~^). The grating was oriented to observe the spectral range A = 2.2194-2.2290 fim in 
a single long-slit spectral order, and a slit position angle of 90° was used. Data were acquired in 
a pair of exposures of 120 s duration each, with the telescope nodded 5" along the slit between 
frames so that object spectra were acquired in all exposures. The BO. 5 star HR 4730 was observed 
for telluric correction of the spectra. Observations of a Th-Ar-Ne hollow-cathode lamp were used 
for wavelength calibrations, and continuum lamp observations were used for flat fields. 



2.5. Data Reduction 



All data were reduced with IRAF. First, object and sky frames were differenced and then 
divided by normalized flat fields. Next, bad pixels were fixed via interpolation, and spectra were 
extracted with the APALL task. Spectra were wavelength calibrated using low-order fits to lines in 
the arc lamp exposures, and spectra at each slit position of each object were co-added. Instrumental 
and atmospheric features were removed by dividing wavelength-calibrated object spectra by spectra 
of early-type stars observed at similar airmass at each slit position. Final spectra were produced 
by combining the spectra of both slit positions for each object and then normalizing them so that 
they had a mean relative flux of 1.0 in each order. 



3. Notes on Individual FUor-Iike Objects 

We now discuss what is already known about the individual FUor-like HHENs, with a focus on 
the properties that indicate episodic variability or other properties that indicate FUor-like behavior. 
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3.1. L1551 IRS 5 



(e.g., iBieging Sz CohenI Il985l ). iRodriguez et al.l (j2003l ) discovered that both binary components 



L1551 IRS 5 was discovered in the near-IR surve y of the Taurus clo ud bv lStrom et al.l (jl976l ). It 
was found to be associated wi th a molecular outflow (ISnell et al.lll98Cl ) and is a close binary system 



drive aligned ionize d bipolar jets from crn wavelength observations. L1551 IRS 5 was confirme d as 
a FUor-like star by lStrom &: StromI (119931 ) after the initial classification by lMundt et al.l (|l985l ). It 
is an IRAS source (04287+1801) with a 12 ij,m flux of ^10 Jy a, nd was later found to be a triple 



system with separations 47 and 13 A U (ILim &: Taka^ 



optical bipolar jet designated HH 154 (jMundt Fried 



kuwa 



19831 ) 



2006 ). It is the driving source of an 



3.2. V883 Ori 



This object was first noted as a faint star illuminating an extensive refiection nebul a (des i gnate d 
IC 430) on photographic plates dating from 1888. In the Ha emission line survey of iHard (|l953l ). 
it was described as being faint with some nebulosity and given the designation H aro 13a. Optical 
spect roscopy of the refiection nebulosity suggested that V883 Ori was a FUor (jStrom Strom 
19931 ). V883 Ori possesses a curving tail of nebulosity, as do many FUors, and is an IRAS source 
(05358-0704) with a 12 /xm fiux of 52 Jy. It has a bolometric luminosity of ~400 Lq and sub- 
mm observations have determined it has an (unresolved) circumstell ar gas+dust mass o f ~0.4 Mq 
(jPent et al.lll998l ). It is thought to be the driving source of HH 183 (jStrom et al.lll986l ). 



3.3. Parsamian 21 



Parsamian 21 was first note d in the catalog of com etary nebulae by iParsamianI (jl965l ). From 
optical images and spectroscopy, IStaude Neckell (jl992l ) concluded that this source was a possible 
FUor and discovere d a small bipolar HH fiow (later labeled as HH 221) eman a ting f rom the star. 
Kospal et al. ( 2007 ) found that the Ha knots first observed by Staude &: Neckei ( 1992 ) were moving 
at velocities 120 - 500 km s~^, and their high resolution near-IR direct and polarimetric images 
reveal a circumstellar envelope, a polar cavity, and an edge-on disk. Parsamian 21 is a ssociated with 
a cold IRAS source (19266+0932) with 12 //m fiux 0.8 Jy and a 100 /xm flux of 15 Jy. IPallv fc Lada 
(|l983l ) found no associated CO outflow from the star. 



3.4. HH 381 IRS 



HH 381 IRS shows prominent optical (jPevine et al. I ll997l ) and near-IR nebulosity (|Connellev et al 



200711 and is an IR AS source (20658+5217) with a 12 ^m flux of 0.3 Jy and 100 fim flux of 11 Jy. 



Devine et all (|l997l ) also considered it the driving source of the HH objects HH 381/382. HH381-IRS 
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was fo u nd to have a near - IR sp ectrum almost identical to that of L1551 IRS 5 bv lReipurth &: Aspin 
(|l997l ). iMagakian et alj (|2007l ) found that the star and associated nebulosity were very faint on 
DSS-1 (1953) plates while bright in DSS-2 (1990) images, and even brighter on recent CCD im- 
ages. Finally, it is not yet conclusive that HH 381 (and in fact HH 380, and 382) originates from 
HH381-IRS since there are several other IRAS sources in the region. 



3.5. HH 354 IRS 



HH354-IRS is lo cated in Lynds 1165 and drives a large-scale HH flow covering ~2.4 pc 
dReipurth et al.lll997l '). [t is associated with an IRAS source (22051-1-5848) and has a weak 12 fim 
flux of 0.3 Jy and a strong 100 fxm flux of 94 Jy. Its bolome tric luminosity is ^120 h g). It was found 
to have a near-IR spectrum very similar to L1551 IRS 5 (jReipurth AspinI 119971) and possess a 
strong sub-mm flux suggesting a gas-|-dust mass of ~ 30 Mq and a CO outflow (jVisser et al.ll2002l ). 



4. Analysis and Results 



Spectra of the classical FUors V1057 Cyg and FU Ori are shown over the 2.203- 2.236 /im and 
the 2.280-2.300 /im ranges in Figured! The GEMINI-South spectrum of FU Ori is not shown here 
because it was superseded by the shown Keck II NIRSPEC spectrum. However, the GEMINI data 
are used in the subsequ ent analysis. Figure [D also sh ows the spectra of the M 4V star GJ 402 (from 
Doppmann et al.ll2005l ). the M2 lab star a Ori f from [Wallace &: Hinkldll996l ). and the veiled (r^ 



1.7) K5-7 Class I protostar p Oph IRS 63 (from Doppmann et al. 20051 ). The spectra of a Ori were 



smoothed so that the CO bandhead has approximately the same slope as the FUors, requiring an 
effective broadening of v sin i ~ 100 km s~^. The a Ori spectra were also artificially veiled by 
= 1.0 so that their features have depths similar to those found in the FUors and FUor-like stars. 
The spectra of the five FUor-like HHENs and a Ori are shown in Figure [2l 



4.1. Spectral Properties 

The M4 V star GJ 402 shows strong neutral atomic absorption lines of Na, Sc, Ti, Fe, 
and Mg species that a re diagnostic of both effective temperature and surface gravity (e.g., see 



Doppmann et al.ll2005l ). The f=2~0 CO bandhead and red- ward vibration-rotation lines are also 
very prominent, and their depths are good indicators of surface gravity when interpreted together 
with the atomic features. The narrow profiles of all lines (except the somewhat gravity-broadened 
Na lines) are indicative of slow, unresolved rotation (v sin i < 17 km s^^). The Class I protostar 
IRS 63 has similar features but they are much weaker due to significant near-IR veiling. They are 



also broad due to the object's considerable rotational velocity, v sin i = 45 km s ^ ( Doppmann et al 



20051 ). The spectra of the M2 lab star a Ori show significant Na, Ti, and CO absorption in ad- 
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dition to many relatively weak features (see [Wallace &: Hinkldll996l ). However, they appear weak 
and broad in Figure [T] due to the artificial veiling and rotational broadening applied. 

The spectra of the two FUors are clearly different from those of the stellar standards and the 
Class I protostar (Figure [I]) . FU Ori shows numerous weak and broad absorption lines throughout 
the 2.280-2.300 fim range whereas the late-type dwarf, supergiant, and protostar show discrete 
features (V1057 Cyg also does not show normal stellar features over its limited spectral coverage). 
The same is true for the spectra of the five FUor-like stars in Figure [21 The v=2-0 CO bandhead 
and red-ward envelope defining that absorption band are certainly the most distinguished spectral 
features in both the FUors and FUor-like stars. 

Optical spectra of FUors have been modeled as being produced in either a circumstellar disk 
(e.g.. lHartmann &: Kenyonlll985l : lKenyon et al.lll988l ) or else in a rapidly rot ating G-type supergiant 
atmosphere surrou nded by a chromosphere and a cooler absorbing shell (e.g., Petrov Sz Herbijll992 : 
Herbig et al.ll2003l ). The near-IR spectrum of FU Ori in the vicinity of the Av = 2 CO bandheads 
and nearby ro-vibration lines has been successfully modeled as arising in a rotating accretion disk 
at radii where the disk photosphe re has an effective temperature and surface gravity characteristic 
of an M-type giant or supergiant ([Kenvon et al.lll988l : iHartmann et al.ll2004l ). Figure [1] shows that 
the larger bandwidth near-IR spectra of FU Ori and V1057 Cyg do share some similarities with the 
artificially broadened and veiled spectra of the M2 supergiant a Ori, but overall they appear to be 
more similar to the spectra of the FUor-like stars in Figure [21 Likewise, the FUor-like stars appear 
to have spectra more similar to FUors than to rotating, veiled late-type stars or Class I protostars. 



4.2. Cross-Correlation Analysis 

We now quantify the similarity of the spectra of the FUors and FUor-like stars with those 
of the stellar dwarfs, stellar giants, and Class I protostars by examining their cross-correlations. 
Figure [3l shows the cross-correlation functions of the 2.203-2.236 fim spectrum of the FUor-like star 
HH 381 IRS against template spectra of FU Ori (Keck spectrum), a Ori (M2 lab), GJ 402 (M4 V), 
and IRS 63 (late-type protostar). Spectra have been shifted slightly to produce symmetrical cross- 
correlation peaks at km s~^, effectively eliminating relative radial velocities. The high central 
peak and symmetrical nature of the cross-correlation in the top panel indicates that the spectra 
of HH 381 IRS and FU Ori have very similar features, but correlations with the giant, dwarf, and 
embedded protostar are poor. 

The broad cross-correlation peaks in all panels of Figure [3l indicate that HH 381 IRS has 
broad spectral features. The peak of the cross-correlation with FU Ori has a half width of ap- 
proximately 100 km s~^, indicating that the 2 //m spectra of these objects are originating from 
regions rotating with velocities v sin i ~ 70 km each (also consistent with the similarly broad 
line widths seen in Figures [H and [2]) . The cross-correlation of HH 381 IRS and the slowly rotat- 
ing M4 dwarf GJ 402 is double-peaked, consistent with the spectrum of HH 381 IRS originating 
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in a rotating disk ([Kenvon et alj Il988l : iHartmann et al.l |2004| ) . This was also seen in the cross- 



correl ation of HH 381 IRS (a nd other FUor-like stars) with the slowly rotating M giant HR 5150 



(from iDoppmann et al.ll2005l . not shown). The cross-correlations with the other stars do not have 
double-peaked structures, likely because of their large rotation velocities (natural in IRS 63 and 
artificially added in a Ori). 

We have computed cross-correlation functions of all FUor-like stars over the 2.203-2.236 /im 
wavelength interval using FU Ori (Keck spectrum), a Ori, HR 5150 (Ml. 5 I II), GJ 402, and the 



protostar IRS 63 as templates. We have computed the iTonrv h DavisI (jl979l ) r-values, a measure 
of the correlation peak height to the average noise, for each correlation function. These values are 
presented in Table [21 with values greater than 3 indicating a significant correlation. 

Table [2] shows that the spectra of the FUor-like HHENs are more strongly correlated with 
FUor spectra than other stars in every case. Also, all FUor-like stars are significantly correlated 
with the spectrum of FU Ori itself (r > 3). The FUor-like stars L1551 IRS 5, Parsamian 21, 
HH 354 IRS, and V883 Ori also show nearly or marginally significant correlations (r ~ 3) with one 
or more other templates, as does FU Ori. This indicates that FUors and FUor-like stars have some 
spectral features in common with giants, dwarfs, or protostars. However, the much more significant 
cross-correlations between FU Ori and FUor-like stars indicate that the spectra of these objects 
are much more similar to each other than to normal stars or embedded protostars. 

In Table [3] we compare the cross-correlation results for each FUor or FUor-like star (over 
2.203-2.236 //m wavelength) with each other in an attempt to measure the similarities of their 
spectral features. Each object is listed in both columns and rows and the cross-correlation values 
are given for each source against each other source. Table entries which cross-correlate a source 
against itself list infinity (oo) as the value. This is not the case, however, for FU Ori, where we have 
cross-correlated the Keck NIRSPEC spectrum taken on UT 2007 Mar 06 (and used in all other 
FU Ori cross-correlations) with the Gemini-S Phoenix spectrum taken on UT 2006 April 03. 

Inspection of the cross-correlation r— values shows that the structure in the spectrum of FU Ori 
correlates extremely well with L1551 IRS 5 (r ~5), V883 Ori (r ~29), Parsamian 21 (r ~22), 
HH 381 IRS (r ~13), and HH 354 IRS (r ~8). This implies that all five sources are spectrally 
very similar to FU Ori. The inter-correlations of these five sources suggests that LI 551 IRS 5 is 
most similar to HH 354 IRS (r ~8), V883 Ori is most similar to FU Ori (r ~29), Parsamian 21 
is most similar to V883 Ori (r ~24), HH 381 IRS is most similar to Parsamian 21 (r ~22), and 
HH 354 IRS is most similar to Parsamian 21 (r ~10). We note also that all cross-correlation 
values are above the r=3 significance threshold indicating that all have significant similarities to 
all others. The significance of the cross-correlation of the two different FU Ori observations is 
very high (r = 9.1) but not infinite. The finite nature of this value is most likely due to the 
limited overlapping spectral range (A = 2.2194-2.2290 /xm) of the 2 observations of FU Ori with 
different instrumentation. Imperfectly corrected instrumental differences may have also contributed 
to reducing the significance of the cross-correlation, and the intrinsic spectrum may have also 
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changed slightly between the two epochs. 



4.3. Near-IR Colors 

The near-IR JHK colors of FUors and HHENs have not previously been compared directly 
and these wavelengths encompass our spectra, so we present a near-IR color-color diagram of 15 of 
these objects in Figure H] (including ones without spectra presented in this paper). The H — K and 
J — H colors in Figure[l]were all computed from 2MASS catalog data, and the plot shows that most 
of the HHENs have colors that are consistent with either reddened classical T Tauri stars (falling 
along the dashed T Tauri locus) or with reddened normal stars with IR excesses; only HH354 IRS 
has colors that are clearly inconsistent with IR excess. FU Ori, several FUors, and a few FUor-like 
stars also have colors consistent with reddened early-type stars (early-type stars are near the origin 
of the plot). However, early- type stars are ruled out and the color degeneracy is broken by the fact 
that the FUors and FUor-like objects generally show strong near-IR CO absorptions, consistent 
with late spectral types typical of T Tauri stars. Therefore most of the FUors and FUor-like 
HHENs have near-IR colors consistent with having circumstellar disks. However, the large range 
of colors makes it difficult to differentiate regular stars, FUors, and T Tauri stars (or protostars) 
by their colors alone. This color range may be produced by some combination of local scattering 
produced by different amounts and distributions of circumstellar material as well as different object 
orientations (i.e., disk / envelope orientations). 



5. Discussion 



The similarity of the FUor and FUor-like spectra and their significant cross-correl ations indi- 

cate th at the spectral lines of these objects must form in physically similar environments. iReipurth &: Aspin 
(|l997l ) had also postulated this based on similarities of low resolution near-IR spectra, and the high 
resolution data presented in this paper confirm that the FUor-like HHENs are more similar to 
FUors than dwarfs, giants, supergiants, or embedded protostars. 

Given the cross-correlation results, the near-IR spectra of FUor-like stars do not appear to 
be produced in normal stellar photospheres. The broad cross-correlation peaks indicate rotational 
velocities v sin i ~ 100 km s~^, niuch larger than ineasured for T Tauri stars or even the most rapidly 
spinning embedded protostars (IDoppmann et al.l l2005l : ICovey et al.l 120061 ). The cross-correlations 
with late- type dwarfs and the embedded protostar IRS 63 are also generally poor (low significance), 
but they are somewhat better with late- type giants and supergiants (Table [2]). This suggests that 
the spectral features of FUor-like stars originate in cool regions with low surface gravities that differ 
from normal stellar photospheres, and the generally good cross-correlations with FUor spectra 
reinforce this. The strong CO absorptions of the FUors and FUor -like stars (Figures [1] ari d [2|) 
indicate relatively low effective temperatures, Teff ~ 3500 K (e.g.. iKleinmann Halll Il986l ). A 
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giant or supergiant star of that effective temperature would have a radius greater than about 40 
Rq and a mass of approximately 6 Mq or more. Such stars have rotational breakup velocities on 
the order of 100 km s~^. Therefore, the CO strengths and projected rotation velocities of FUor and 
FUor-like spectra (from cross-correlation and CO bandhead widths; Fig. [2] and [3|) are consistent 
with giants rotating at breakup. However, the marginal cross-correlation significance with giants 
or supergiants (Table [2]) indicates that the FUors and FUor-like stars are not likely to be ordinary 
giants or supergiants. 

The combined visible and IR spectra o f FUors are reasonably well fit by models of line formation 
in active circumstellar accretion disks (e.g. jKenvon et al.lll988l : lHartmann et al.ll2004l ): the decrease 
in their line widths (and implied rotational velocities) at longer wavelengths is recognized as perhaps 
the best evidence of line formation in disks. The double-peaked nature of the cross-correlation of 
FUors and FUor-like near-IR spectra with those of late type dwarfs (Figure [3]) is also consistent 
with a disk rotational velocity profile, and the significant near-IR excesses exhibited by most of 
these objects (e.g., §4.3 and Fig. H]) are also consistent with hot circumstellar disk material. 

If the spectra of FUors and FUor-like HHENs are indeed produced in accretion disks, then 
why are they not more similar to the spectra of embedded protostars with active disks? Even 
Class I protostars with large veil ing and relatively high accretion rates M ~ 1 0~^Mf7^ yr~^ do not 
have FUor-like spectra (see also I Greene &: Ladal l2002l : iDoppmann et al.l l2005l ) . In a circumstellar 
disk model atmosphere, high accretion rates are needed to induce a temperature gradient with 



a hot midplane and coole r exterior to produce a spectrum similar to a FUor (e.g., iKenyon et al 



19881 ). ICalvet et al.l (|l99ll ) found that an irradiated accretion disk model with mass accretion rate 
of 1.6 X 1O~^M0 yr~^ reproduced the near-IR spectrum of FU Ori well, including its CO and 
H2O absorptions. This accretion rate is approxir nately 2 orders of ma gnitude higher than that 
of embedded low mass protostars such as IRS 63. iGreene fc Ladal (j2002l ) estimated that the mass 
accretion rate of the embedded protostars YLW 15A was 2.3 x 10~^Mq yr~^, and the same modeling 
would produce a very similar (within ~ 50%) value for IRS 63 since the two protostars have virt ually 
identical effective temperatures, near-IR veilings, rotational velocities (IDoppmann et al.ll2005l ) and 
IRAS fluxes. It appears that higher accretion rates are needed to produce luminous disks with 
absorption features similar to those found in FUors. Future discoveries and observations of weaker 
FUors with lower accretion rates and more active protostars with higher accretion rates would 
be valuable in constraining further the accretion rate at which a circumstellar disk develops a 
strong enough temperature gradient to produce absorption lines that dominate those of the stellar 
photosphere. 

The very high degree of similarity of the FUor-like and FUor spectra (and their good cross- 
correlations) suggest that FUor activity may be more common than indicated from the small number 
of known classical FUors alone; more young stars have spectra likely to originate in luniinou s 



accretion disks. Taking only the list of 20 classical and FUor-like objects by Abraham et al 



and adding the two other HH sources found to have FUor-like spectra by iReipurth &: AspinI ( 



(2004) 



19971 ) 



(HH 381 IRS and HH 354 IRS; also confirmed in this work) yields a sample of 22 objects. Of 
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these, five are widely acknowledged to be bona fide classical FUors due to their observed outbursts 
and spectral properties (FU Ori, V1057 Cyg, V1515 Cyg, V1735 Cyg, and V346 Nor). Therefore 
FUor-like stars may outnumber FUors. 

However, FUors are certainly not ubiquitous. Very few FUors or stars with FUor-like spectra 
are found in regions of clustered star formation such as Tau-Aur, p Oph, NGC 1333, IC 346, 
or the Orion Nebula Cluster. Instead, they are mostly found in regions of low star formation 
activity. This implies that either FUors are much less frequent or have many fewer outburst cycles 
in high density regions. One possible explanation is that interactions between stars in clusters could 
disrupt disks, perhaps eliminating much FUor activity during the Class or Class I protostellar 
evolu tionary phases. Alternat ively, FUor events may be triggered by the presence of a companion 
star (jBonnell &: BastienI Il992l ) , in which case a dense cluster environment might accelerate the 
orbital evolution, causing FUor eruptions in clusters to occur mostly during the Class or Class I 
stages (iReipurth &: Aspinll2004l ). However, these are certainly not definitive explanations, and the 
causes of FUor events must be understood better before the inhomogeneity in the distribution of 
FUors and FUor-like stars can be understood. 
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Fig. 1. — Ncar-IR spectra of two FUors, a Class I protostar, and late- type stellar standards. The 
FUors do not strongly resemble the other spectra except for having strong CO and some absorptions 
near the Na lines. 
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Fig. 2. — Spectra of five FUor-like stars and the M2 lab star a Ori. The spectra of a Ori have 
been smoothed so that the CO bandhead has approximately the same slope of the FUor-like stars, 
V sin i ~ 100 km s~^. 
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Fig. 3. — Cross-correlation functions of the 2.203-2.236 //m spectrum of the FUor-hke star 
HH 381 IRS. The high central peak and symmetrical nature of the cross-correlation in the top 
panel indicates that the spectra of HH 381 IRS and FU Ori have similar features, but correla- 
tions with the giant, dwarf, and embedded protostar are poor. Spectra have been shifted slightly 
in wavelength (and therefore radial velocity) to produce symmetrical cross-correlation peaks at 
km s~^. 
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Fig. 4. — Near-IR color-color plot of classical FUors (filled star symb ols) and FUor-like s tars 
(open star symbols). Several stars posited to be FUor-like stars (e.g., see lAbraham et al.ll2004l . for 
references) but not otherwise analyzed in this work are included. J UK magnitudes were obtained 
from the 2MASS database, and uncertainties in colors are generally smaller than the symbol size. 
Most FUors and FUor-like stars have near-IR colors that are consistent with reddened T Tauri 



stars with disks (classical T Tauri locus of iMeyer et al .1119971 . shown as the dashed line) but not 
reddened late- type dwarfs (open circles and their locus) or giants (filled dots and their locus). The 
dot-dashed line indicates reddening equivalent to = 10 mag. 
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Table 1. Journal of Observations 



Source 


a(2000) 


5(2000) 


Observed Date 


Int. Time 


S/N 


Observatory 




(hh mm ss.s) 


(° ' ") 


(UT) 


(minutes) 






FUor-like HHENs: 












L1551 IRS5 


04 31 34.1 


+18 08 05 


2001 Nov 6 


20.0 


90 


Keck 


V883 Ori 


05 38 18.1 


-07 02 27 


2007 Mar 6 


1.0 


140 


Keck 


Parsamian 21 


19 29 00.7 


+09 38 39 


2001 Jul 7 


33.0 


220 


Keck 


HH 381 IRS 


20 58 21.4 


+52 29 27 


2001 Jul 7 


21.0 


190 


Keck 


HH 354 IRS 


22 06 50.7 


+59 02 49 


2001 Jul 7 


21.0 


60 


Keck 


FUors: 














FU Ori 


05 45 22.4 


+09 04 12 


2007 Mar 6 


1.0 


120 


Keck 








2006 Apr 3^ 


4.0 


160 


Gemini- S 


V1057 Cyg 


20 58 53.7 


+44 15 29 


1999 Aug 30"^ 


13.0 


300 


IRTF 








1996 Sep 2^ 


2.0 


250 


IRTF 



^Single order spectrum of the 2.2194-2.2290 /xm region with spectral resolution R ~ 40,000. 
This spectrum is not shown in the figures of this paper (as it is largely superseded by the Keck 
spectra), but it is used in the cross-correlation analysis. 

^Spectrum of the 2.29353 fim v=0-2 CO bandhead region with spectral bandwidth ~ 57A. 

^Spectrum of the 2.2075 fim Na line region with spectral bandwidth is ^ 55A. These spectra 
were originally published in Greene &; Lada (1997). 
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Table 2. Cross-correlation r- values 



Object 


FU Ori 


a Ori 
M2 lab 


HR 5150 
M1.5 III 


GJ 402 
M4 V 


IRS 63 
protostar 


L1551 IRS 5^ 


4.1 


2.5 


3.2 


2.3 


3.5 


V883 Ori 


28.7 


4.3 


4.0 


3.8 


3.8 


Parsamian 21 


21.6 


3.1 


3.2 


2.9 


3.1 


HH 381 IRS 


13.3 


1.6 


2.5 


1.7 


2.3 


HH 354 IRS 


7.5 


3.0 


3.0 


1.9 


3.3 


FU Ori 


oo 


3.4 


2.9 


3.3 


3.1 



'^Thc 2.223 H2 emission line was artificially removed from the spectrum of L1551 IRS 5 before 
computing its cross-correlations. 



Table 3. Cross-correlation Between Objects (r-valucs) 



Object 


L1551 IRS 5^ 


V883 Ori 


Parsamian 21 


HH 381 IRS 


HH 354 IRS 


FU Ori 


L1551 IRS 5^ 


00 


5.0 


6.0 


6.6 


8.1 


5.1 


V883 Ori 


5.0 


00 


25.2 


18.8 


8.7 


28.7 


Parsamian 21 


6.0 


25.2 


00 


21.5 


9.6 


21.6 


HH 381 IRS 


6.6 


18.8 


21.5 


00 


6.7 


13.3 


HH 354 IRS 


8.1 


8.7 


9.6 


6.7 


00 


7.5 


FU Ori 


5.1 


28.7 


21.6 


13.3 


7.5 


9.1^ 



'^The 2.223 //m H2 emission line was artificially removed from the spectrum of L1551 IRS 5 before 
computing its cross-correlations. 

^The 2007 Mar 6 Keck NIRSPEC spectrum of FU Ori, shown in Figure 1 and used for all FU Ori 
cross-correlations, was cross-correlated with a spectrum of FU Ori acquired on 2006 April 03 with 
the GEMINI-S telescope (See Table 1). 



